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> ■ 1 INTRODUCTION 



ABSTRACT 

Evolution characteristics of a Kerr black hole (BH) are investigated by considering 
coexistence of disc accretion with the Blandford-Znajek process (the BZ process) and 
magnetic coupling of the BH with the surrounding disc (MC process) . (i) The rate of 
extracting energy from the rotating BH in the BZ process and that in MC process 
are expressed by a unified formula, which is derived by using an improved equivalent 
circuit, (ii) The mapping relation between the angular coordinate on the BH horizon 
and the radial coordinate on the disc is given in the context of general relativity 
and conservation of magnetic flux, (iii) The power and torque in the BZ process are 
compared with those in MC process in detail, (iv) Evolution characteristics of the BH 
and energy extracting efficiency are discussed by using the characteristics functions of 
BH evolution in the corresponding parameter space, (v) Power dissipation on the BH 
horizon and BH entropy increase are discussed by considering the coexistence of the 
above energy mechanisms. 
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It is well known that magnetized accretion disc of black hole 
(BH) is an effective model in astrophysics. The interaction 
between a Kerr BH and the surrounding magnetic field has 
been used not only to explain high energy radiation and 
jet production from quasars and active galactic nuclei, but 
also as a possible central engine for gamma-ray bursts (Rees 
1984; Frank, King & Raine 1992; Lee, Wijers & Brown 2000, 
hereafter LWB). Confined by the magnetic field in the inner 
region of the disc, the magnetic field lines frozen previously 
in the disc plasma will deposit on the horizon in company 
with the accretion onto the BH, and the magnetized disc 
becomes a good environment for supporting the magnetic 
field on the horizon. Blandford and Znajek (1977) proposed 
firstly that the rotating energy and the angular momentum 
of a BH can be extracted by the surrounding magnetic field, 
and this energy mechanism has been referred to as the BZ 
process, in which the BH horizon and the remote astrophys- 
ical load are connected by the open magnetic field lines, and 
energy and angular momentum are extracted from the ro- 
tating BH and transported to the remote load. 

There still remain some open problems on the BZ pro- 
cess, and one of them is how to estimate the ratio of the 
angular velocity of the magnetic field lines to that of the 
BH horizon. Macdonald and Thorne (1982 hereafter MT82) 
argued in a speculative way that the ratio will be regu- 



lated to about 0.5 by the BZ process itself, which corre- 
sponds to the optimal value of the extracting power with 
the impedance matching. However, as argued by Punsly and 
Coroniti (1990), it is hard to understand how the load can 
conspire with the BH to have the same resistance and satisfy 
the matching condition, since the load is so far from the BH 
that it cannot be casually connected. 

Recently some authors (Blandford 1999; Li 2000a, 
2000b; Li & Paczynski 2000) argued that with the existence 
of the closed field lines a fast rotating BH will exert a torque 
on the disc to transfer energy and angular momentum from 
the BH to the disc. Henceforce this energy mechanism is 
referred to as magnetic coupling (MC) process. Compared 
with the BZ process, the load in MC process is the surround- 
ing disc, which is much better understood than the remote 
load though the magnetohydrodynamics (MHD) of the disc 
is still very complicated. 

Some works have been done to discuss the influence of 
the BZ process on the evolution of BH accretion disc (Park 
& Vishniac 1988; Moderski & Sikora 1996; Lu et al. 1996; 
Wang et al. 1998, hereafter WLY). However MC effects in- 
volving the closed field lines were not taken into account in 
the previous works. In fact the closed field lines should ex- 
ist on the horizon as well as the open field lines, and disc 
accretion should coexist with the BZ and MC process. In 
this paper evolution characteristics of a Kerr BH are inves- 
tigated by considering coexistence of disc accretion with the 
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BZ and MC process (henceforce DABZMC). 

In order to facilitate the discussion of the evolution of 
the central BH surrounded by the magnetized accretion disc 
we make the following assumptions: 

(i) The disc is perfectly conducting and the magnetic 
field lines are frozen in the disc. The magnetosphere is sta- 
tionary, axisymmetric and force-free outside the BH and the 
disc (MT82); 

(ii) The magnetic field is so weak that its influence on 
the dynamics of the particles in the disc is negligible, and 
the BH has an external geometry of Kerr metric; 

(iii) The disc is thin and Keplerian, lies in the equato- 
rial plane of the BH with the inner boundary being at the 
marginally stable orbit. 

This paper is organized as follows. In Sec. II the rate of 
extracting energy from the rotating BH in the BZ process 
(henceforce the BZ power) and that in MC process (hence- 
force MC power) are expressed by a unified formula, which 
is derived by using an improved equivalent circuit based on 
MT82. Our result for the BZ power is consistent with that 
derived in LWB. In Sec. Ill a mapping relation between the 
angular coordinate on the BH horizon and the radial coordi- 
nate on the disc is given in the context of general relativity 
and conservation of magnetic flux. In Sec. IV The power and 
torque in the BZ process are compared with those in MC 
process in detail. In Sec.V evolution characteristics of the 
BH and energy extracting efficiency are discussed by using 
the characteristics functions in the corresponding parame- 
ter space. In Sec. VI we discuss the power dissipation on the 
BH horizon and the BH entropy increase, and show that the 
excess rate of change of BH entropy does come from the to- 
tal power dissipation on the BH horizon in the BZ and MC 
process. Finally, in Sec. VII, we summarize our main results. 



2 DERIVATION OF BZ POWER AND MC 
POWER 

Two mechanisms of magnetic extraction of energy from a ro- 
tating BH are involved in our model, i.e., the BZ process and 
MC process. The configuration of the poloidal magnetic field 
of the BH magnetosphere for the coexistence of DABZMC 
is shown in Fig.l, where the angular coordinate of the open 
field lines is assumed to vary from to 9m, and that of the 
closed field lines from 9m to ir/2. 9m is a parameter indicat- 
ing the angular boundary between the open and the closed 
field lines, and the BZ and MC effects vanish as 9m = and 
9m = 7r/2 , respectively. 

Considering that both mechanisms arise from the rota- 
tion of the Kerr BH relative to the surrounding magnetic 
field, we expect that the BZ power and MC power can be 
expressed by a unified formula. In order to do this we pro- 
posed a model, in which the remote load in the BZ process 
and the load disc in MC process are incorporated into a 
unified load with resistance Zl and angular velocity £7l . 
Based on MT82 an improved equivalent circuit is shown in 
Fig. 2, in which a series of loops correspond to the adjacent 
magnetic surfaces consisting of the field lines connecting the 
BH horizon and the load. 

In the i-loop of Fig.2 segments PS (characterized by the 
flux \P) and QR (characterized by the flux *P + A*) repre- 
sent two adjacent magnetic surfaces, and segments PQ and 




Figure 1. The configuration of the poloidal magnetic field of the 
BH magnetosphere. 
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Figure 2. An improved equivalent circuit for a unified model for 
the BZ and MC process. 



RS represent the BH horizon and the load sandwiched by 
the two surfaces. AZ H and AZ L (hereafter subscript "i" is 
omitted) are the corresponding resistances of the BH hori- 
zon and the load, respectively. Aeh = (A 1 !/ /2tv)Qh and 
Ael = — (Aty /2-k)Q,l are electromotive forces due to the 
rotation of the BH and the load, respectively. The minus 
sign in the expression of Ael arises from the direction of 
the flux. Q.h in the expression of Aeh is the angular veloc- 
ity of the BH horizon and reads (G=c=l): 



Q.h = a,/(2r H ), r H = M (1 + q) , q = y/l - a 2 



(1) 



where a* is the BH spin, which is related to BH mass M and 
angular momentum J by a* = J/M 2 , and m is the radius 
of the BH horizon. By using the BH magnetosphere model 
given in MT82 we obtain the following equations: 



f£ aE-dl = AV H = (A*/2tt) (luh - SIf) 
= (A*/2tt) {Q h - Of) = IAZh 

J*a(E + v L xB)- dl = AV L + (A*/2tt) (u>l 
= (A*/2tt) (fi F - fl L ) = IAZ L 



AVr = 



aE-dl= (A*/2tt) (fi F - u) L ) 



(2) 



(3) 



(4) 



where AVh and AVl are the potential drops measured by 
"zero-angular-momentum observers" (ZAMOs) on the BH 
horizon and on the load, respectively. ZAMOs is a family 
of fiducial observers defined by Bardeen et al (1972). vl is 
the load's velocity relative to ZAMOs. uih (— » fifj) and lol 
are the ZAMO angular velocities on the BH horizon and the 
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load, respectively. Of is the angular velocity of the magnetic 
field lines, a = (dr/dt) ZAMO is the "lapse function", which 
is written as follows: 



a = pVA/£ (5) 
where E, p and A are the Kerr metric parameters and read 



/ 2 . 2\ 2 2 a • 2 Q 2 2 2 2 

(r +a ) —a Asm 8, p =r +a cos 



A = r 2 + a 2 - 2Mr 



(6) 



Incorporating equations (2) and (3) we obtain 

Ae H + Ae L = / aE-dl + a [E + v l x B) ■ dl 

J P J Ft 

= I (AZ H + AZ L ) (7) 

Dividing equation (2) by equation (3), we derive the relation 
relating the three angular velocities, £Ih, and Of as 
follows: 

Q F = (QhAZl + Q l AZh)/(AZ h + AZ L ) (8) 

Assuming Ql = in the BZ process, we have 

Q F = VL H AZ L /{AZ H + AZ L ) (9) 

Considering that the load disc consists of plasma of perfect 
conductivity and the magnetic field lines are frozen in the 
disc plasma, we have AZl = and 



Q F = = ^1 



M( X 3 + a t ) 



(10) 



where fl D is the angular velocity of the disc a t the place 
where the magnetic flux penetrates, and \ = \/r/M is di- 
mensionless radial coordinate of the disc. The current I in 
each loop can be expressed by 

/ = A j H + A l L =(A<g/27Q 



AZ H + AZ L 



AZ H + AZ L 



(A*/2tt) ■ 



AZ L 



(11) 



where equation (8) is used in deriving the right-hand side 
(RHS) of equation (11). The current / on the BH horizon 
feels Ampere's force, and the torque exerted on the annular 
ring of the width Al — pA6 is 



AT = ujBhIAI = (A*/2tt) 2 (Qh - Q, F )/AZ H 



(12) 



where A* = Bh^it-cdAI, Bh is the magnetic field on the 
horizon, and vj = (S/p) sin# is the cylindrical radius of the 
horizon. Therefore the extracting power between the two 
adjacent magnetic surfaces is 

ap = n F AT = (a*/2tt) 2 n F (n H - n F )/Az H (13) 

where 

AZ H = RhAI/(2ttvj) = 2 P A6/vj (14) 

and Rh — 4tt/c — 377ohm is the surface resistivity of the 
BH horizon. Considering that the BH horizon consists of 
two hemispheres and integrating equations (13) and (12) 
over the angular coordinate from 9 = to 8m we obtain the 
total BZ power and torque as follows: 



Pbz/Po = 2al 



' k (1 - Jfc) sin 3 OdO 
2- (l-q) sin 2 <9 



(15) 



Tbz/T q = 4a» (1 + q) 



f 0M (l-k); 
Jo 2-(l- 



sin d 9d9 



q) sin 2 9 



(16) 



where k = £I f /Q,h and 



Po = (B 2 H ) M 2 w B 2 A M'i x 6.59 x 10 4 

(b%) M 3 w BlM'i x 3.26 x 10 4 ' g ■ cm 1 ■ s" 



/ ?- A/f3 r?2 ]{/f3 ^ y o oc 1 n47 „ „_2 „ — 2 V V 



In equation (17) (i?lr) is the average value of B% over the 
BH horizon, B4 and Ms are yj (Bfj) and M in the units 
of 10 4 gauss and 10 s Mq, respectively. Hereafter {B^ is 
regarded as a constant on the horizon, and the sign of the 
average value is omitted. Taking 9m = tt/2 and k = 0.5 in 
equation (15) we obtain the optimal BZ power as follows: 

P° p * imal /P = Q- 1 (arctanQ - a,/2) , 

Q EE x J(l-q)/(l + q) (18) 

This result is exactly the same as derived by Lee et al, who 
pointed out that the BZ power had been underestimated 
by a factor ten in previous work (LWB). The corresponding 
optimal BZ torque can be expressed by 



T o P timal/p o = 4Q -2 ( arctan Q _ at / 2 ) 



(19) 



Similarly, taking £l F — and integrating equations (13) 
and (12) over the angular coordinate from 9m to tv/2, we 
obtain total MC power and torque as follow: 



Pmc/Po = 2a* 



7I-/2 



13(1-13) sin 3 9d9 
2 -(l-q) sin 2 9 



Tmc/T = 4a, (1 + , 



w/2 



(1-/3) sin 3 9d9 
2 -(l-q) sin 2 9 



where 



/3ee n D /n„ = 



2(l + g) 
a» (x 3 + a*) 



(20) 



(21) 



(22) 



Inspecting equations (15), (16), (20) and (21), we find the 
expressions for Pbz and Tbz are almost the same as those 
for Pmc and Tmc except that k is replaced by (3. From 
equation (9) we know that the parameter k = Q F /Q H is 
uncertain due to lack of the knowledge about "the remote 
astrophysical load" , while (3 can be determined by equation 
(22) as a function of a* and X- 



3 MAPPING RELATION BETWEEN BH 
HORIZON AND MC REGION OF DISC 

In order to calculate Pmc and Tmc we should first deter- 
mine the mapping relation between the BH horizon and the 
disc. Considering the flux tube consisting of two adjacent 
magnetic surfaces "i " and "i + 1 " as shown in Fig.l, we 
have A^h = A^d required by continuum of magnetic flux, 
i.e., 



BH2n(vop) r=rH d9 = —Bd2tt (vjp/\TK) g _ 7T ^ 2 dr 



(23) 



where Bo is the poloidal component of the magnetic field 
on the disc and 



(tup) r =r H = (£sin#) r = rjf = 2MrHsin# 
(wp/VA) e=n/2 = S/\/A = a _1 p 



(24) 
(25) 
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where 



«=(pVa/e) = 



2 X 



(26) 



1 + a^XmtC- 2 + 2a2 x - 6 4_ 3 

Following Blandford (1976) we assume that Bo varies as 

B D oc C™ (27) 

where £ = r/r ms is a dimensionless radial parameter defined 
in terms of the radius r ms of the marginally stable orbit 
(Bardeen et al 1972). 

Assuming that the inner boundary of the MC region is 
located at the inner edge of the disc, and considering the 
balance of the magneic pressures between the horizon and 
the inner edge of the disc (Ghosh & Abramowicz 1997), we 
assume 



Bh = Bd 



(28) 



Incorporating equations (27) and (28) we have 

B D = BhC" (29) 

Substituting equations (24), (25) and (29) into equation (23) 
we have 



sinddd = — G (a*;£, n) d£ 
where 



(30) 



G(a*;£,n) 



XmsS 



2(1 + ?) 



1 + aix™Z- 2 + 2^X^~ 3 
l-2 X m4- 1 +a*Xm 4 4" 2 



(31) 



Integrating equation (30) we have 
ens Dm — cos (i = — I G (a*; £, n)d£ 



cos 9 m 



I 

Jtout 

. at 9 -- 
G(a*;£,n) 



Setting £ = £ in = 1 at 6 — tt/2, we have 



(32) 



(33) 



Thus the mapping relation between the angular coordinate 
on the horizon and the radial coordinate on the disc is de- 
rived in the context of general relativity and conservation of 
magnetic flux as follows: 

re 

cos 6= I G(a,;£,n)d£ (34) 



-I 



It is noted that equation (33) provides a constraint to the 
outer boundary of the MC region, which is represented by 
£ out = r out /r ms , provided that the parameters a*, 9m and 
n are given. £ out turns out to be a function of n, varying 
monotonically for the given a* and 9m as shown in Fig.3, and 
it behaves as a function of a* varying non-monotonically for 
the given n and 9m as shown in Fig. 4. Since £ out is limited, 
the parameter n should be less than one critical value n c 
corresponding to infinite £ ou t- The curves of n c versus a* for 
the given values of 9m are shown in Fig. 5. 



4 COMPARISON OF POWER AND TORQUE 
IN THE BZ PROCESS AND MC PROCESS 

In order to compare the relative strength of MC power to the 
BZ power and torque we define the ratio of Pmc to p°P^ tmal 



s 





Figure 3. The curves of t; ut versus n for the given values of 8m 
with 6m = (solid line), it/ 6 (dashed line) and it/ A (dotted line); 
(a) a* = , (b) a, = 0.998. 




Figure 4. The curves of t; ut versus a* for 8m = and < a» < 
0.998 with n = 1.1 (solid line), n = 1.5 (dashed line), and n = 3 
(dotted line). 
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Table 1. The three special values of 7 and the values of the 
concerning parameters 



n 


7 = 


= -1 


7 


= 


7 = 


- ^max 


1.1 


aZ 1 = 


0.2583 


al 2 = 


0.3286 


aj m 


= 0.7997, 




af = 








7max 


= 0.9278 


1.5 


0.2581 


aT = 


0.3283 


a I 


= 0.8002, 












Tin ax 


= 0.9287 


3.0 


71 
a* = 


0.2572 


a* = 


0.3269 


a I 


= 0.8021, 












7max 


= 0.9325 



Table 2. The three special values of r and the values of the 
concerning parameters 



Figure 5. The curves of n c versus a* for the given values of 6m 
with 8m = (solid line), n/6 (dashed line) and it /4 (dotted line). 



-1 




Figure 6. The curves of 7 (solid line) and r (dashed line) versus 
the BH spin a* for n = 1.1 and 0.2 < a* < 0.998. 

and that of T MC to T"""" 1 " 1 as follows, 



= Pmc/P\ 



optimal 
BZ 



lalQ 



(arctanQ — a*/2) 



r/2 



0(1 -/3)sin 3 6d6 
2 - (1 - q) sin 2 6 



^_ rp ^ J rpoptimal 

alQ 



(arctanQ — a*/2) 



r/2 



(1-/3) sin 3 6d6 
2- (1-g) sin 2 6» 



(35) 



(36) 



Setting 9m = in equations (35) and (36), we have the 
curves of 7 and r versus a* for the given parameter n, and 
those versus n for the given a* as shown in Fig. 6 and Fig. 7, 
respectively. 

The values of a* corresponding to the three special val- 
ues of 7 (7 = —I, and ^ ma x) and to those of r(r = —1, 
and T max ) are calculated for the given values of n as shown 
in Table 1 and Table 2, respectively. 

From the above calculations on 7 and r we obtain the 
following results: 

(i) The signs of Pmc and Tmc are the same as the 
corresponding 7 and r, and they change from the negative to 
the positive at a~l 2 and a T t 2 respectively. It is shown in Table 



1.1 


o-l 1 


= 0.1988 


< 2 


= 0.3278 


a rm 


= 0.8045, 












7~max 


= 0.7452 


1.5 


o-l 1 


= 0.1987 


< 2 


= 0.3275 


< m 


= 0.8056, 












7"max 


= 0.7481 


3.0 




= 0.1980 


< 2 


= 0.3260 


a rm 


= 0.8118, 












7"max 


= 0.7614 



1 and Table 2 that a» 2 is a little less than aj 2 . Therefore both 
energy and angular momentum are transferred from the BH 
to the disc as a» > a* 2 , and transferred from the disc to the 
BH as a» < aj 2 . It is interesting to notice that the transfer 
direction of energy is opposite to that of angular momentum 
for a T2 < a, < a" 2 , i.e., energy is transferred from the disc 
to the BH, while angular momentum from the BH to the 
disc in this value range of the BH spin. 

(ii) Both 7 and r vary non-monotonically as a« for the 
given parameter n, and attain their maxima at a» m and 
a™, respectively; 

(iii) The maxima of Pmc and Tmc are all less than 
the corresponding r BZ and 1 BZ , respectively, while 
the absolute values of the former two are greater than the 
corresponding values of the latter two as a* is less than 0Z 1 
and al , respectively. 

(iv) As shown in Fig. 7, both 7 and r increase very little 
as the increasing parameter n, which implies that both Pmc 




Figure 7. The curves of 7 (solid line) and r (dashed line) versus 
the power-law index n for a* = 0.5 and 1 < n < 6. 
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and Tmc are insensitive to the power-law index n of the 
radial coordinate £. 



5 EVOLUTION CHARACTERISTICS AND 
ENERGY EXTRACTING EFFICIENCY 

Based on conservation of energy and angular momentum the 
basic evolution equations of the Kerr BH in the coexistence 
of DABZMC are written as follows: 

dM/dt = E ms M D - P BZ - Pmc (37) 

dJ/dt = L ms M D -T BZ -Tmc (38) 

Incorporating equations (37) and (38), we have the evolution 
equation for the BH spin as follows: 

da,/dt = M~ 2 (L ms M D - T BZ - T MC ) 

-1M- x a*{E ms M D -PBz-PMc) (39) 

where Md is the accretion rate of rest mass, E ms and L ms 
are specific energy and specific angular momentum corre- 
sponding to the inner edge radius r ms , respectively. 

Since the magnetic field on the BH is supported by the 
surrounding disc, there are some relations between Bh and 
Md- As a matter of fact these relations might be rather 
complicated, and would be very different in different situ- 
ations. One of them is given to investigate the correlation 
between BH spin and dichotomy of quasars by considering 
the balance between the pressure of the magnetic field on 
the horizon and the ram pressure of the innermost parts of 
an accretion flow (Moderski, Sikora & Lasota 1997), i.e., 

— Pram ^ pC ~ M D /(4irr 2 H ) (40) 

From equation (40) we assume the relation as 

M D = B 2 H M 2 (1 + q) 2 /2 (41) 

However equation (41) is derived without MC effects. By 
considering the transfer of the energy and angular momen- 
tum between the BH and the disc equations (41) is modified 
as follows: 

M D = B 2 H M 2 (1 + q) 2 /2 + [1-8 (9 M - tt/2)] (M D ) mc (42) 

where [Md ) is the MC correction to the accretion rate 
at the inner edge of the disc, and it reads 

4B 2 H M 2 Q( x 3 ms +a,) 2 (l-f3) I ' l + ajx^^al^E ,.oN 
(«k s +4a.) V l-2 x -^ x -i ^ 

The derivation of equation (43) is given in Appendix A. To 
exclude (Md) from the case without magnetic coupling 

V / mc 

we define 8 {9 m — k/2) as follows: 

/Hi oTo:%2 <*> 

Thus the equation (42) is applicable to the coexistence of 
DABZMC. 

Setting Md = in equation (42), we can derive a value 
of the BH spin, at op ~ 0.4677, at which the accretion will 
stop at the inner edge of the disc. Substituting equation (42) 
into equations (37)— (39), we have 



dM/dt = / (a, , Uj ) M D (45) 

dJ/dt = Mh(a*,Uj) M D (46) 

da- t /dt = M' 1 g(a,,u j )M D (47) 
where 

g(a*,Uj) — h(a t ,Uj) — 2a„ / (a, , Uj ) (48) 



In the above equations Uj = 8M,k,n is used for simplicity 
to represent the three parameters other than a*. f(a,,Uj) 
and g(a*,Uj) are referred to as the characteristics functions 
(CFs) of BH evolution. 

5.1 Evolution characteristics of Kerr BH 

First we discuss the evolution characteristics of the Kerr BH 
in the coexistence of disc accretion with the BZ process for 
6m = tt/2. From equations (37) — (39) we have the CFs as 



follows: 

/ (a» , k) = E ms -P BZ /M d (49) 

h(a t ,k)=M- 1 (L ms -T BZ /M D ) (50) 

g(a*,k)=h(a*,k)-2a,f(a*,k) (51) 



From equations (45) and (47) we find that the signs of 
dM/dt and da*/dt are the same as those of f(a,,k) and 
g (a*, k), respectively. As shown in Fig.8 we have the curve 
represented by g (a, ,k) — in the 2-dimension space 
consisting of the parameters a* and k, which divides the 
parameter space into two regions. And each black dot 
with arrowhead is referred to as a representative point 
(RP), which represents one BH evolution state. We can use 
RP's displacement in the parameter space to describe BH 
evolution visually. According to the positions of RPs we 
have two evolution modes of the BH, and the details are 
given as follows. 

(I) Region I: /(a,, k) > and g{a,,k) < 0. 

Mode I: RPs in Region I always move towards the 
left until they arrive at the boundary curve g(a t ,k) = 0, 
i.e., BH will reach its equilibrium spin a% q in long enough 
evolution time. RPs in region I represent the BHs with 
increasing mass and decreasing spin. 

(II) Region II: f(a*,k) > and g{a,,k) > 0. 

Mode II: RPs in Region II always move towards the 
right until they arrive at the boundary curve g (a*, k) = 0, 
i.e., BH will reach its equilibrium spin at q in long enough 
evolution time. RPs in region II represent the BHs with 
increasing mass and increasing spin. 

Next we shall discuss BH evolution in the coexistence 



of disc accretion with MC process. Substituting 9m = 
and equation (42) into equations (37) — (39), we have the 
corresponding CFs as follows: 

/ (a» , n) = E ms - Pmc / M d (52) 
h(a„n) = M' 1 (L ms - T MC / M D ) (53) 
g (a„ ,n) = h (a« , n) — 2a*/ (a„ , n) (54) 



We can also describe BH evolution by using the displace- 
ment of RPs in 2-dimension parameter space consisting of 
the parameters a* and n as shown in Fig. 9, which is divided 
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Figure 8. Parameter space for BH evolution in the BZ process 
with < k < 1 and 0.7 < a* < 1. 
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Figure 9. Parameter space for BH evolution in MC process with 
1 < n < 6 and 0.4598 < a» < 0.4606. 



into two regions by the curve <; (a*, n) = 0. In this case the 
description of the BH evolution characteristics and the two 
evolution modes are almost the same as that given in the 
BZ process except that f(a*,k) and g(a,,k) are replaced 
by / (o», n) and g (a*, n), respectively. Inspecting Fig. 8 and 
Fig. 9, we have the following results: 

(i) The CFs in the corresponding parameter space can 
be used to describe the evolution characteristics of the BH 
in a visual way. If we know the initial position of RP in the 
parameter space, we can determine the evolution character- 
istics of the BH immediately. 

(ii) The effects of the parameter k and n on the BH 
evolution in the BZ and MC process can be easily found in 
the corresponding parameter space. From Fig. 8 we find that 
af increases as the increasing k, which is consistent with 
our previous work (WLY). On the other hand, we find from 
Fig. 9 that at q increases very little as the decreasing n. Our 
calculations show that the values of a eq are very near to and 
a little less than af ° v ^ 0.4677. 

(iii) As shown in Fig. 5, the value range of n is limited 
and confined by conservation of magnetic flux. Compared 
with the parameter k, the effects of n on the variation of at q 
are very weak for the confined value range of n. 



5.2 Energy extracting efficiency of BH accretion 
disc 

Energy extracting efficiency of BH accretion disc can be in- 
vestigated also by using the CFs. From equations (45) and 
(47) we have 



dM/M — [f (a,,iij)/g(a,,Uj)]da, 



(55) 



As is well known, a Kerr BH can be characterized by its 
mass M and spin a,. If a Kerr BH evolves from the initial 
state (a,to,Mo) to the state (a,,M), the BH mass can be 
expressed by 



M = M exp 



g{a*,Uj) 



da„ 



(56) 



The efficiency corresponding to the evolution process of BH 
can be defined as 



_ M D -(M- Mq) _ Mj Mp - 1 
Vp ~ W D ~ ~ M D /Mo 



(57) 



which is the efficiency of converting accreted mass into ra- 
diation energy during the evolution process, and Mo is the 
accreted rest mass and expressed as a function of a, by using 
equations (47) and (56) as follows: 



Mr 



Mo 



■ exp 



a, f ( a»,u 



,0 9 



(a'„Uj) 



da, 



da*(58) 



The efficiency corresponding to each evolving state of the 
BH is defined as 
dM D - dM 



dM D 



(dM/dt)/M D = 1- f{a*,Uj) (59) 



From equation (37) we find that efficiency rj a is equal to the 
sum of the following three terms: 



Vs 



Vda + Vbz + Vmc 



(60) 



where ijda = 1 — E ms , r) B z = Pbz/M d and -que = 
Pmc/Md are the contributions of disc accretion, the BZ 
process and MC process, respectively. From equation (47) 
we know that the equilibrium spin at q is determined by 



g(a*,Uj) = 

Incorporating equations (55) — (59), we have 
V P \ a ,^= Mm (1-A/B) 

In equation (62) we have 

ca. / (a' t ,uA i 



(61) 
(62) 



and 



ra, J \ a„,Uj\ ( 

A = exp / — r— ; r-da* 

A.o g(a,,Uj) 

pa, 1 pa, f I a,, Uj I 

B= — -exp / — - — T da, 

Ja, 9( a *i U i) Ja, g{a*,Uj) 

By using L' Hospital limit rule for oo/oo type we have 

^P\a,^af = 1 -/( ffl * 9 > U j) =Vt q {<C,Uj) (63) 

where r;f 9 represents the efficiency rj a at at q . Equation (63) 
implies that the efficiency r\ v is exactly equal to rf q , provided 
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Figure 10. The curves of a% q and rf s q varying as the concerning 
parameters, (a), (d) 6m = 7r/4, < k < 1, n = 1.1; (b), (c) 
e M = tt/4, k = 0.5, 1 < n < 6; (c), (f) < d M < tt/2, n = 1.1, 
A: = 0.5(solid line), k = 0.788(dashed line); k = 0.96(dotted line). 



that a BH evolves from any spin a*o to its equilibrium spin 
a q . Incorporating equations (42), (59), (61) and (37) — (39), 
we obtain the curves of a% q and varying as the concerning 
parameters are shown in Fig. 10. 

Comparing Fig. 10a, b, c with Fig.lOd, e, f we find the 
following results: 

(i) The curves of a% q and rf q varying as the concerning 
parameters look very alike except that rf q varies as k non- 
monotonically with its maximum (r/s q )max « 30.5% at a% q w 
0.4611 for k w 0.8998 as shown in Fig.lOd. 

(ii) Both at q and rf q decrease monotonically as the in- 
creasing n for 9m = 7r/4 as shown in Fig. 10b and Fig.lOe. 
Considering that the outer boundary £ out increase mono- 
tonically as the increasing n, we infer that the more concen- 
trated is the magnetic field on the central region of the disc, 
the wider is the MC region on the disc, the stronger is the 
braking effects of MC process on the BH spin, and the lower 
are the at q and r)l q . 

(iii) As shown in Figs. 10c and lOf the variations of a eq 
and rf q as 9m display different characteristics for different 
values of parameter k: They might increase as 9m (dotted 
line) or decreases as 9m (solid line) or almost remain un- 
changed with 9m (dashed line). It is not difficult to under- 
stand these results by considering that 9m is the angular 



boundary between the open and closed field lines on the BH 
horizon: The different variations of at q and rf q as 9m come 
from the different contributions due to the BZ and MC pro- 
cess. For example, af and rf q will increase as 9m , if the 
contributions of the open field lines in unit angular coor- 
dinate are greater than those corresponding to the closed 
ones, while they will decrease as 9m in case that the former 
is less than the latter. 

In summary we can calculate the efficiency in an evolu- 
tion process and that corresponding to each evolving state 
of the Kerr BH by using the equations (57) and (59), re- 
spectively, and both of which are expressed in terms of the 
CFs. 



6 BH ENTROPY CHANGE IN COEXISTENCE 
OF DABZMC 

It was shown that the excess rate of change of entropy of the 
central BH in the BZ process comes from the total power 
dissipation on the BH horizon (WLY). It is attractive for us 
to extend this result to the coexistence of DABZMC. As is 
well known, the BH temperature Th and entropy Sh can be 
expressed as (Thorne, Price & Macdonald 1986, hereafter 
TPM) 

q 



T H = 



,Sh = 2nM'(l + q) 



(64) 



47rM(l + i 

Incorporating the basic evolution equations (37), (38) and 
(39), we have the following equations, 



(65) 



T H dS H /dt = dM/dt - Q H dJ/dt = (E ms - 
Q, H L ms )M D + (n H T B z - Pbz) + (VhTmc - Pmc) 

equation (65) is exactly the mathematical formulation of the 
first law of thermodynamics for a Kerr BH (TPM). Obvi- 
ously, the first term on RHS of equation (65) is the con- 
tribution due to disc accretion. It is easy to prove that the 
second and the third term arise from power dissipation on 
the BH horizon due to the BZ and MC process, respectively. 
By using the improved equivalent circuit as shown in Fig. 2 
and equation (11) we have 



IA £l 



I 2 AZ h + I 2 AZ L 



Incorporating equations 
- (A*/27r) Q l we obtain 



- IAsl (66) 
(9), (11), (13) and Ae L = 



AP = (A*/2tt) 



2 Qf (Off 



Of) 



AZ, 



I l AZ L - IAe L 



(67) 



From equations (66) and (67) we find that IAeh consists 
of two parts: the extracting power AP and the dissipated 
power APdsp- The extracting power is transferred to the 
load including both the remote astrophysical load and the 
load disc. For the BZ process Ael = 0, we have 

, 2 Q, F (Q H -n F ) 



AZ h 



(68) 



2 SId (fiff 



fir.) 



AP = AP BZ = I AZl = (A*/2tt 

For MC process AZl = 0, we have 

AP = APmc = -I Ae L = (A*/2tt) 

Incorporating equations (11), (12) and (13), we can express 
the dissipated power APbsp on the BH horizon as 

APdsp = I 2 AZ H = (A*/2tt) 2 {Q h - Vt F f / AZ H 



(69) 
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Q H AT - AP 



(70) 



Integrating equation (70) over the angular coordinate 9, we 
obtain the total dissipated power on the horizon, which is 
exactly equal to the sum of the second and the third term 
on RHS of equation (65), i.e., 



Pdsp = (Pdsp) bz + (Pdsp) mc 



(71) 



where (Pdsp) bz = [QhTbz — Pbz) and (Pdsp) mc = 
(FIhTmc — Pmc) are the total power dissipation due to the 
BZ and MC process, respectively. Combining equations (15), 
(16), (20) and (21) we have 



(Pdsp), 



(Pdsp) mc 



QhTbz — Pbz 



2aiB 2 H M 2 



f 

Jo 



(1 - fcfsin 3 9d9 
2-(l-q) sin 2 9 



QhTmc — Pmc 

f-T/2 



2a 2 B 2 H M 2 



r 



(1- /3) 2 sin 3 6d9 
2 - (1 - q) sin 2 9 



(72) 



(73) 



Incorporating equations (65), (72) and (73), we express the 
rates of change of BH entropy corresponding to disc accre- 
tion, the BZ and MC process as follows: 

(dS H /dt) DA = Tx 1 (E ms - n H L ms ) Mo 



(dSu /dt) BZ = Tf, 1 (QhTbz — Pbz) 
(dSu /dt) MC — T H L (CIhTmc — Pmc) 



(74) 
(75) 
(76) 



From the above analysis we conclude that the excess rate of 
change of the BH entropy does come from the total power 
dissipation on the BH horizon in the BZ and MC process. 

The ratios of the rates of BH entropy change can be 
written by 

(dSn/dt)BZ 



Rbzda (a* , 6m , k) — 

QhTbz — Pbz 
(E ms — n_f/L ms ) Mo 



(dS H /dt)i 



Pmc da (a* , 9 m ,n) = 



(dS H /dt) MC 
(dS H /dt) DA 



MC 



i Ml 



(77) 



(78) 



where the four parameters k, a,, 9m and n are involved. 
Letting one of them vary and the rest be fixed, we have the 
curves of Rbzda and Rmcda versus a„ and 6m as shown 
in Fig. 11, and the following results are obtained: 

(i) As shown in Figs. 11a and lib Rbzda increases 
monotonically as a* for the given values of k and 9m, while 
Rmcda varies non-monotonically and attains a minimum 
approach zero at a* w 0.3305 for the given values of n and 
9m- From the discussion in Sec. IV we know that the min- 
imum of Rmcda arises from the neighboring a* 2 and a r2 
corresponding to Pmc = and Tmc = 0, respectively. 

(ii) As shown in Figs. 11a and lib we have Rmcda > 
Rbzda for < a* < 0.2694 and 0.4373 < a* < a% q , which 
implies that the contribution to Pdsp and BH entropy in- 
crease due to MC process dominates over that due to the 
BZ process in the above value range of the BH spin. 

(iii) As shown in Fig. 11c the ratio Rbzda increases 
monotonically as 9m for the given values of a, and k, while 




0.05 



(a) 




0.45 



(b) 




Figure 11. The curves of Rbzda (solid line) and 
Rmcd A (dashed line) versus the concerning parameters, (a), (b) 
< a* < 1, 6 M = Tr/ 4 : k = 0-5, n = 1.1; (c) < 6 M < tt/2, 
a„ = 0.4, k = 0.5, n = 1.1. 



Rmcda varies non-monotonically as 6m for the above given 
values of a, and n. Rmcda is greater than Rbzda for < 
9 M < 0.6782, while the former is much less than the latter 
for 0.6782 < 9m < tt/2. It is obvious that the sum of the two 
is much less than unity, which implies that the contribution 
to BH entropy increase is dominated by disc accretion in the 
coexistence of DABZMC. 
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7 SUMMARY 

In this paper the evolution characteristics of a Kerr BH 
and some related issues, such as energy extracting efficiency 
and entropy change on the BH horizon, are investigated by 
considering coexistence of DABZMC. By using an improved 
equivalent circuit a unified expression for the BZ power and 
MC power are derived. Starting from the conservation laws 
of energy and angular momentum, we obtain the basic evo- 
lution equations of the BH and the CFs in terms of four 
parameters related to our model, i.e., a*, 9m, k and n. We 
find that the evolution characteristics of the Kerr BH can 
be well described by using the CHs and the RPs in the cor- 
responding parameter space. 

It turns out that the main difference between the BZ 
process and MC process lies in the two kinds of loads: the 
remote load with an unknown parameter k for the BZ pro- 
cess, and the load disc with the parameter f3 for MC process. 
Parameter k is uncertain, while f3 is thoroughly determined 
by the BH spin and the place where the magnetic flux pen- 
etrates. The effects of MC process on the evolution of the 
BH and the efficiency of BH accretion disc can be discussed 
in virtue of the parameters o*, 9m and n, into which the 
parameter j3 is merged. 

In this model the mapping relation between the BH 
horizon and the disc is derived by assuming a power-law of 
the magnetic field varying as the radial coordinate of the 
disc. It is shown that the power-law index n is related to 
the outer boundary parameter £ ou t- And the MC correction 
to the accretion rate is considered by using the conservation 
law of angular momentum. Only the accretion rate at the 
inner edge of the disc being involved in the evolution equa- 
tions of the BH, the correction to the accretion rate is made 
by taking the corresponding values at r ms in our simplified 
model. 

In summary this model provides an analytic approach to 
the evolution characteristics and related physical quantities 
of a Kerr black hole surrounded by magnetized accretion 
disc. 
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APPENDIX A: DERIVATION OF 
EQUATION (43) 

In the case without magnetic coupling between the BH and 
the disc, accretion is produced by the internal viscous torque 
T V i 3 of the disc. By using the conservation law of angular 
momentum the accretion rate Md is related to T„i s by the 
following equation (Frank, King & Raine 1992): 



M D d(r 2 n D )/dr = dT m8 /dr 



(Al) 



where dT v i s /dr is the contribution due to internal viscous 
torque, always transporting angular momentum outward in 



the disc. Taking MC effects into account, we think equation 
(Al) should be modified by 

- M D d(r 2 n D )/dr = dT vls /dr - dT MC /dr (A2) 

where dTMc/dr is the contribution due to MC torque, and 
the minus sign arises from the consideration that the con- 
tribution of dT v i a /dr should be counteracted by dTMc/dr 
in case Qd < Qh- Comparing equations (Al) and (A2), 
we express (Mr)) mc , the contribution of MC effects to the 
accretion rate, as follows: 



on 



dT MC /dr 



d{r 2 Q, D ) d(r 2 D. D )/dr 



(A3) 



Substituting equations (6), (10), (12), (14) and (34) into 
equation (A3), we have 



(Md), 



4B 2 H M 3 a„(l + q)(l - j3) sin 2 9 



2 - (1 - q) sin 2 
2( X 3 W 2 



sX 2 (x 3 +4a,; 



G(a*;£,n) 



(A4) 



From equation (A4) we find that the MC correction de- 
pends on the place on the disc where the MC torque acts. 
Substituting the boundary values such as £ = £i n = 1 
(i.e.,r = r ms — Mxms) and 9 — ty/2 into equation (A4) 
we obtain equation (43), which is the MC correction to the 
accretion rate at the inner edge of the disc. In derivation the 
factor 2 is given for the two faces of the disc, and Bjj is also 
taken as the average value over the BH horizon. 
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